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ABSTRACT 

The solar photosphere is a partiaUy ionized medium with coUisions between electrons, 
various metallic ions and neutral hydrogen playing an important role in the momentum 
and energy transport in the medium. Furthermore, the number of neutral hydrogen 
atom could be as large as 10^ times the number of plasma particles in the lower pho- 
tosphere. The non-ideal MHD effects, namely Ohm, Ambipolar and Hall diffusion can 
play an important role in the photosphere. We demonstrate that Hall is an important 
non-ideal MHD effect in the solar photosphere and show that Hall effect can signif- 
icantly affect the excitation and propagation of the waves in the medium. We also 
demonstrate that the non-ideal Hall dominated inhomogeneous medium can become 
parametrically unstable, and it could have important ramification for the photosphere 
and chromosphere of the sun. The analysis hints at the possibility of solar photosphere 
becoming parametrically unstable against the linear fluctuations. 

Key words: Sun: Photosphere, MHD, waves. 



1 INTRODUCTION 

The solar atmosphere is a partially ionized medium with 
fractional ionization {X^ — n^jrin) - a ratio of the electron 
(rie) to the neutral (n„) number densities, changing with 
the altitude. In fact, the lower solar atmosphere - the pho- 
tosphere 500 km) is a weakly ionized medium where col- 
lisions between electrons, metallic ions and neutrals are re- 
sponsible for the dynamical behaviour of the medium. With 
increasing altitude, the fractional ionization, which is 10~^ 
in the photosphere drops to 10~^ in the chromosphere at 
1000 km. Even at an altitude ~ (2000 - 2500) km, the neu- 
tral number density dominates the plasma density by an 
order of magnitude (Vernazza et al. 1981). Thus, the solar 
atmosphere is a partially ionised mixture of plasma and neu- 
tral gas with very weakly ionized gas in the photosphere and 
fully ionized plasma in the upper chromosphere and beyond. 

Clearly, a partially ionized mixture of gas does not be- 
have like fully ionized plasma and hence, ideal magnetohy- 
drodynamics (MHD) description of the dynamics of solar 
photosphere at best is a very crude approximation. The in- 
vestigation of Alfvenwave in the ideal MHD has been very 
popular topic in the solar physics owing to its appeal to 
the solar coronal heating (lonson 1978; Priest 1987; Goed- 
bloed & Poedts 2004). These waves can be easily excited 
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in the medium due to easy availability of vast reservoir of 
energy such as the convective gas motions. Only a tiny frac- 
tion of this energy, if carried by the wave to higher altitude, 
will suffice to heat the corona to high temperatures (Parhi 
et al. 1997). Implicit in this picture is the assumption that 
non-ideal MHD effects are unimportant in the photosphere. 
However, collision effects may sever ly affect the propagation 
of the wave in the solar atmosphere (De Pontieu & Haeren- 
del 1998; Goodman 1998, 2004; Leake et al. 2005; Vranjes 
et al. 2007; Arber et al. 2007; Vranjes et al. 2008). 

The importance of non-ideal MHD effect is related to 
the question of how well the magnetic field is coupled to the 
neutral matter. These effects can be quantified in terms of 
plasma-Hall parameter /3j 



ft 



(1) 



a ratio of the j^'^ particle cyclotron frequency tJcj = eB /rrij c 
(where e ,B , raj , c denots the electron charge, magnetic 
field, mass and speed of the light respectively) to the sum of 
the plasma-plasma, and plasma - neutral, j^jnCoUision fre- 
quencies. For electrons — Ven + v^i + Vee ~ 2 i^en + and 
for ions Vi = Vin -\- va- It is clear that whenever /3j ^ 1, the 
j^^ plasma particle will be strongly tied to the ambient mag- 
netic field of the medium whereas in the opposite case, par- 
ticles will not feel the Lorentz force and thus, can be treated 
as weakly magnetized or unmagnetized. For example, when 
ions and neutrals are strongly coupled, the relative drift of 
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the "frozen-in" ions (/3i ^ 1) against the sea of neutrals 
may cause the diffusion of the magnetic flux - the so called 
ambipolar diffusion. When the ions and neutrals are mov- 
ing together (/3i < 1), and there is relative drift between the 
electrons and ions (i.e. electrons arc "frozcn-in" the field but 
ions are not), the Hall diffusion ^ becomes important. When 
even electrons are not "frozen-in" the field, i.e. /3e ^ 1, the 
Ohmic dissipation becomes important. Therefore, the inves- 
tigation of the non-ideal magnctohydrodynarnic (MHD) ef- 
fect in the solar atmosphere should be investigated in differ- 
ent parameter windows: (I) ambipolar - when the magnetic 
field can be regarded as frozen in the plasma and drifts with 
it through the neutrals (II) Ohmic when neutrals stops 
the ionized particle from drifting with the field, and (III) 
Hall - wheu electrons are well coupled or partially coupled 
to the field and ions are partially or completely decoupled 
from the field. Since Hall parameter (3j will vary with the 
altitude, we shall anticipate that solar atmosphere will be in 
different non-ideal MHD regimes at different altitudes. 

As noted above, non-ideal MHD effects are dependent 
upon the ambient physical parameters of the medium and 
their relative importance, and consequently, their role in the 
excitation and propagation of the waves can be gauged by 
the plasma-Hall parameter. In table 1, we give some rele- 
vant parameters for the solar photosphere (Vernazza et al. 
1981) and corresponding plasma-Hall parameter. We have 
assumed that all ions are hydrogen ions although metal- 
lic ions generally dominate the cold photosphere. However, 
numbers for metal ions are speculative as the cross-section 
is not known. Assuming cross section a = 30 x ain for the 
metallic ions {rrii = 30 mp, here mp = 1.67 10"^* g is pro- 
ton mass) gives similar Hall beta parameters as for the ions 
(table 1). Therefore, we shall refrain from giving Hall pa- 
rameter separately for the metallic ions. As is clear from 
the values of fUj at different altitude, the Ohmic dissipation 
will be dominant process at the surface of the photosphere. 
With the increasing altitude. Hall diffusion becomes the im- 
portant mechanism for the flux distribution in the medium. 
Although numbers are not given for 250km, Hall staxts to 
competes with the Ohmic diffusion at this altitude. For a 
10^ G field. Hall will dominate the Ohmic diffusion at the 
surface. The ambipolar diffusion is unimportant in the lower 
photosphere (^ 1500 km) if B = lOG, since A < 1. Only 
at higher altitude, typically between (1500 — 2500) km, am- 
bipolar diffusion becomes important although Hall is only 
half as small as ambipolar in this region. Since neutral rmm- 
ber density plummets rapidly beyond 2500 km, the role of 
ambipolar diffusion diminishes in the upper chromosphere. 
However, since cause of the Hall effect is due to the sym- 
metry breaking between electrons and ions with respect to 
the magnetic field (Pandey & Wardle 2008), Hall will con- 
tinue operating at this height as well albeit on a smaller 
scale. This conclusion appears counterintuitive since Hall 
MHD in a fully ionized plasma generally operates when the 
frequencies of interest is larger than the ion-cyclotron fre- 
quency over ion-inertial scale 5i = VAi/i^d, a ratio of the ion- 
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Alfvon speed {vai ~ B/^An pi, pi = rniUi is the ion mass 
density) to the ion-cyclotron frequency. For a 10 G mag- 
netic field, the cyclotron frequency is 10® , 10^ and lO'^ ,s~^ 
at h = 0,515 and 1065 km respectively. Further the Hall 
scale, in a fully ionized plasma. Si « 2, 44 and 71 cm, corre- 
sponding to to /i = , 515 , and 1065 km respectively, is very 
small. Clearly, the spatial and temporal scales of Hall MHD 
derived from a fully ionized, two component electron ion 
plasma suggest that the Hall effect is unimportant. There- 
fore, fully ionized Hall MHD description of solar atmospheric 
plasma is irrelevant. Understandably, the role of Hall MHD 
of a fully ionized plasma in the solar atmosphere has been 
overlooked by the solar community. 

The importance of Hall diffusion in angular momentum 
transport in protoplanctary discs (Wardle 1999; Balbus & 
Terquem 2001; Sano & Stone 2002a,b; Salmeron & War- 
dle 2003, 2005; Pandey & Wardle 2006) indicates that Hall 
operates on different spatial and temporal scales in a par- 
tially ionized plasmas. Recently, Pandey & Wardle (2008) 
developed a single fluid formulation for the partially ionized 
plasma and showed that the spatio-temporal scale over which 
Hall operates in a partially ionized plasma is very different 
from a fully-ionized Hall-MHD. Defining Hall frequency as 
(Pandey & Wardle 2008) 

UlH = — OJci , (2) 
P 

where p — pn + pi = rn„ n„ -\- rrii Ui is the mass density of 
the bulk fluid with m„ , nii as the mass and n„ , Ui as the 
number density of the neutral and ion fluids respectively, 
we can demand that if dynamical frequency of the system is 
larger than the Hall frequency, i.e. ujh ^ w, then Hall effect 
will be important. Furthermore, Hall in a partially ionized 
medium will not operate on the ion-inertial scale as in a fully 
ionized plasma, but over (Pandey & Wardle 2008) 

L.^(l + ^) 5i. (3) 

Assuming m; = mn, we see from Table 1 that Hall frequency 
varies between 10 — 10® Hz and thus for fluctuations occur- 
ring at higher than loh , Hall will play an important role over 
Lh ~ few hundred meters to few kms in the medium. We 
note that unlike high frequency Alfvenwave which damps 
in the solar medium (Leake et al. 2005), the excitation of 
low-frequency normal modes due to collisional coupling will 
propagate un-damped in the medium. 

We investigate the propagation of waves in the solar 
photosphere in the Hall regime using recently developed gen- 
eral set of equations by Pandey & Wardle (2008) applicable 
to the partially ionized plasmas. The effect of dissipative dif- 
fusion on the wave properties has been investigated in the 
past (Tanenbaum & Mintzer 1962; Kulsrud & Pearce 1969). 
The wave dissipation in such a medium is dependent not 
only on the ion-ncutral collision frequency but also on the 
fractional ionization of the medium (Kumar & Roberts 2003; 
Pandey & Wardle 2008). The wave damping have recently 
been studied in the context of spicules dynamics (Pontieu 
& Haerendel 1998). In the lower solar photosphere however, 
since Hall will dominate all other diffusive processes the re- 
sultant low frequency ion-cyclotron and collisional whistler 
wave will be the normal mode of the medium (Pandey & 
Wardle 2008). 
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Table 1. The neutral mass density p = rUn nn, the ratio of ion to neutral mass density, the ion and electron 
collision frequencies and plasma Hall parameters at different heights of the solar atmosphere are given in 
the table. The ion mass mj = rtij, and the magnetic field B = 10 G has been assumed. 



h(km) 


rin (cm 3) 


ne/n„ 




(Hz) 




(Hz) 




(Hz) 




(Hz) 


ft 


/9e 





1.2 • IQ^'' 


5.3 • 10-4 


1.2 


•lO^ 


2.2 


•10^ 


3.6 • 


ioi° 


1.3 


•109 


10-^ 


10-2 


515 


2.1 • 10^5 


1.2 • 10-4 


1.8 


•10^ 


2.1 


■ 10^ 


5.2 


•108 


1.3 


•10^ 


10-3 


10° 


1065 


1.7 ■ 10^3 


10-2 


1.7 


• 10^ 


5.1 


• 104 


5.1 


• 10« 


3.1 


• 10^ 


10-1 


10 


1515 


1012 


6-10-2 


1.1 


• 104 


3.3 


• 104 


3.2 


• 10^ 


2.0 


• 10^ 


10° 


102 


2050 


7.7 ■ 10i° 


5 • 10-1 


8.8 


• 10^ 


1.6 


• 104 


2.6 


• 104 


9.6 


• 10^ 


10° 


102 


2543 


10^ 


1.2 • 10° 


0.9 


• 102 


1.1 


• 10° 


2.6 


• 10^ 


6.6 


•lOi 


10 


10^ 



We investigate the linear and nonlinear wave properties 
of the medium in the present work. We show that since such 
a medium is inherently dispersive, the balance between the 
dispersion and nonlinearity leads to DNLS equation. The 
basic set of equations and the linearized dispersion relation 
is discussed in Sec. II . In Sec. Ill we first discuss the para- 
metric instability and then describe the nonlinear equation. 
In Sec. IV discussion of the results and a brief summary is 
presented. 



2 BASIC MODEL 

The solar photosphere is a weakly ionized medium consisting 
of electrons, ions, neutrals. The single fluid description of 
such a plasma has been given in the past (Cowling 1957; 
Braginskii 1965). We shall use the single fluid description 
given by Pandey & Wardle (2008). 

The continuity equation for the bulk fluid is given as 

dp 
dt 

where p = pe + pi + pn ~ Pi + Pn is the bulk fluid density 
and V is the bulk velocity, v = {piVi + pnVn)/p with p,, n 
and p„, Vn as the mass density and bulk velocities of the ion 
and neutral fluids respectively. The momentum equation is 

given as 



+ V-(pv) = 0, 



(4) 



dv „ „ J X B 

p— = -VP+ . 

at c 



(5) 



Here J = rie (vi — Ve) is the current density, B is the mag- 
netic field and P = Pe + Pi + Pn is the total pressure. The 
induction equation is 



dB 

at 



Vx 



„N 4 TT ?7 4:71 riH ^ „ 

{v X B) J — J X B 



c 

4iTr]A 



(j X B^ 



X B 



(6) 



where B — B/B, and the Ohmic (77), ambipolar (tja) and 
Hall {rin) diffusivity are 

c2 D^B^ 



V = 



4: TV Pi Vi, 



,m = 



cB 



(7) 



Here D = Pn/p- We sec from above Eq. (6) that the ra- 
tio of the HaU (H) and the Ohm (O) gives H/O ~ /3e 
and the ratio between ambipolar (A) and Hall (H) gives 
A/H ~ Pi. In the photosphere, Z? ^ 1, /3i < 1 and 
thus ambipolar diffusion can be neglected. In Fig.l, we 
show the ratio oi A/H and H/O for the collision frequen- 
cies from table 1 for B = 100 G field. The figure suggest 
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Figure 1. The ratio of Hall to Ohm (/3e), and ambipolar 
to Hall (£>2 /J.) is given for a B = 10^ G field. 



that between < /i ^ 620 km in the solar atmosphere. 
Ohm will dominate both non dissipative Hall and dissipa- 
tive ambipolar diffusion; between 620 km ^ ft 1160 km 
Hall dominates Ohm as well as ambipolar whereas between 
1160 km ^ /i ^ 2500 km ambipolar will dominate Hall. 
Clearly these estimates are based on the constant magnetic 
field strength. In Fig. 2, we show that ambipolar to Hall 
ratio changes with the changing field strength. For exam- 
ple, if B = 10 G, the Hall will dominate the ambipolar in 
1000 km ^ ft ^ 1600 km and only after 1600 km ambipolar 
becomes more important than Hall. Below 100 km Ohmic 
dissipation dominates all other form of diffusion. With in- 
creasing magnetic field strength, the doiiiaiu of Ohmic dis- 
sipation shrinks and Hall starts operating at much closer 
to the surface. Therefore, we may say that for sufficiently 
strong field B > IkC, Hall diffusion will operate in the 
photosphere whereas ambipolar difi^usion will operate in the 
chromosphere. Since Hall scale is dependent upon the frac- 
tional ionization of the medium (Pandey & Wardle 2008), 
with increasing altitude, the Hall scale Lh = Vh/va will 
shrink and finally approach the ion inertial scale for a fully 
ionized plasma. 

The role of coUisional effects in the photosphere and 
chromosphere has been investigated in the past (Goodman 
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Figure 2. The ratio of ambipolar and Hall term is shown 
in this figure for varying magnetic field strength. 



2004) using the conductivity tensor approach. We note that 
in the conductivity tensor approach plasma inertia is ig- 
nored while deriving generalized Ohms law, and MHD waves 
(which require non-zero plasma inertia!) are studied using 
such a generalized Ohms law. Clearly to overcome this logi- 
cal inconsistency, a general set of equations Eqs. (4), (5), and 
(6) have been derived by Pandey & Wardle (2008). These set 
of equations can bo closed by prescribing a thermodynamic 
relation between pressure and mass density. In order to com- 
pare the results of present formulation with the previous 
investigations e.g. (Goodman 2004), we note that the gener- 
alized Ohms law J = a\\ E^^ + a± E± +aH E x B/B (here 
E = E + vxB is the field in the bulk frame, E\\ = E ■ B, 
Ej_ = B X (e X B^ where B = B/B) can be inverted to 
yield 



E' = ± 



(8) 



Here a± = ^/a'j, + a'jj and J± = B x{J xB). The parallel 
(c7||), Hall (an) and Pedersen (ap) conductivities are given 
by Wardle & Ng (1999) 



(Tp 



c e Ui 
B 



B 



ceiii 
B 



(/3e + I3i) 



f3e 



Pi 



1+/ 



1 + 



(9) 



(10) 



(11) 



While writing above expressions for the conductivity ten- 
sors, we have assumed that plasma is quasineutral, i.e. 
ne ~ rii. Taking curl of Eq. (8) and using cVxE = —dtB 
we get the induction equation (6) except that 77 , rjH and 
r?A will be expressed as a combination of ff\\,ap and au- 
Assuming /3e ^ 1 and /3i ~ 1, we note that 



where a = Ba j [cen). In the /3e S> 1 and /?, ~ 1 limit , 
the ratio of ambipolar and Hall terms from Eq. (8) can be 

written as 



(13) 



which is same as derived from induction Eq. (6) except for 
the factor. This difference arises because in the conduc- 
tivity tensor approach, ion inertia is set to zero and thus 
D = pnj {pi + pn) = 1 since pi = 0. Clearly the induc- 
tion equation of Pandey & Wardle (2008) is more general 
than the conductivity tensor approach which is valid only 
for weakly ionized plasma and is ill suited to describe the 
transition regions where plasma inertia is important. The 
ratio of Hall to Ohm from Eq. (8) becomes 



— 

_2 Pe J 



(14) 



which is the same as derived from cq. (6). Clearly in the 
weakly ionized limit Pandey & Wardle (2008) and conduc- 
tivity approach gives identical results. It has been suggested 
in the past that Pedersen current dissipation could play an 
important role in the chromosphere (Goodman 2004). We 
see from Fig. 2 that ambipolar and Hall effect dominates the 
lower chromosphere (particularly for a strong magnetic field) 
and thus, the dissipative Pedersen current heating could be 
a plausible mechanism for coronal heating. 

Eqs. (4), (5), and (6) together with an equation of state 
P = ci p forms a closed set. We note that the coUisional 
dissipation will cause the loss of energy in the lower solar 
atmosphere and therefore, generally a proper energy equa- 
tion should be used for a more realistic modelling of the 
physical processes (Arber et al. 2007). However, in order to 
keep the description simple, we shall use P = p to inves- 
tigate the wave properties of the medium. In order to study 
waves in the photosphere, we shall assume a homogeneous 
background with no flow (vq — 0). The equations (4), (5) and 
(6) after linearizing around Bo, Po Po and Fourier analysing 
exp {iujt — ik ■ x) becomes 



u6p — pk-Sv = 0, 

-1 



uidv ■■ 



4tt p 



{k- B)dB - [ — ) {SB ■ B) k 



(15) 
(16) 



Here lu'^ = lo^ — k'^ c^. Defining k ■ B — cos 6, along with 
Alfvon frequency loa = kvA, for SB we get 
2 

(w^ - uj\ cos^ e)dB = ^ uil (^SB -B^ (^B-k cos 6^ 

—IK r]H oj cos e (fe X (5b) . (17) 

Only Hall and convoctivo terms have been retained while 
linearizing the induction Eq. (6). After some straightforward 
algebra, following dispersion relation can be derived from 
equation (17) 



(w^ — ui\ cos^ 6) X 



■ uiA cos^ ^1 = k'^ ri%iJ^ cos^ . (18) 



; , ap w Pi , an ^ 1 , 



(12) 



Since compressional mode will not affect the transverse 
mode, we shall drop ci term in the above dispersion re- 
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lation. Then writing rin = lo'a/'^h, we can recast above 
Eq. (18) as 



2 2 2fl,/WX2 a 

01=0;^ COS a ± ( ) loa cos a . 



(19) 



This dispersion relation acquires a famihar form (cf. Wardle 
& Ng 1999, Eq. 25) when wave is propagating along the 
ambient magnetic field (0 = 0) 

In the lja ^ i^h limit, i.e. when the dressed ions (with 
effective mass m* — p/rii) gyrates faster than the time over 
which magnetic fluctuation propagates, the dispersion rela- 
tion (18) gives the familiar Alfven wave 



2 2 a 

; uja cos tf . 



(20) 



We note that these waves propagate in the neutral medium. 
The propagation of these low frequency Alfven waves in a 
neutral medium occurs due to coUisional coupling of the 
plasma particles with neutrals. The electrons, ions and neu- 
trals are well coupled by the collisions and partake in the 
oscillations together. 

The propagation of waves in a collision dominated 
medium appears counter intuitive since one would expect 
that collision will dissipate the energy and thus damp- 
ing of the fluctuations will occur. This picture is valid 
only when wave frequencies are comparable or greater than 
the effective neutral-ion collision frequency, i.e. kvA ^ 
\/pn/pi (/3e/l + D (3e) Vni implying that the neutral is not 
hit often enough by the ions to participate in the oscil- 
lations (Pandey & Wardle 2008). Therefore, only fluctu- 
ations of certain wavelength will disappear due to coUi- 
sional dissipation, and waves of wavelengths A exceeding 
'^cutoff ~ DvA/vni Can propagate in the medium. 

Since magnetic restoring force due to field deformation of 
wavelength A acts only on the plasma particles, the neutrals 
at a distance A apart can respond simultaneously to this 
restoring force only if the communication time between the 
neutrals (~ \/va) is smaller than the effective collision time 
tc ~ Dv~l. Only when tc exceeds ~ \/va, the wave in the 
medium will damp. Using table 1, we see that that since 
^cutoff ^® ^-^^ ^™ at h = and 7 m at ft = 1000 km, the low 
frequency Alfven wave of largo wavelengths will propagate 
in the partially ionized solar photosphere without damp- 
ing. Therefore we can say that the solar photosphere sup- 
ports the excitation and propagation of the low frequency 
Alfven waves in the predominantly neutral medium where 
the inertia of the fluid is carried by the neutral fluid and the 
magnetic deformation is felt by the plasma particles. 

When loh tOA, i.e. when the dressed ion gyration 
period is slower than the propagation time of the magnetic 
fluctuations in the medium, the dispersion relation (18) can 
be analysed in low cu <C i^a and high lua <S lv frequency 
limits. In the low frequency limit, we get 

2 2 2 a 

U! ~ LOff cos . 



In the high frequency loa ^ OJ limit. 



(^VcOS^^. 



(21) 



(22) 



The mode described by Eq. (21) is the modified electrostatic 
ion cyclotron mode. Since w ~ ujh ~ 10 - lO'' Hz, solar 
photosphere can support ion cyclotron waves, except when 
the direction of wave propagation is almost transverse to 



the ambient magnetic field. Equation (22) is the dispersion 
relation for the whistler mode. In the lower photosphere, 
the whistler is high frequency branch whereas ion-cyclotron 
frequency is the low frequency branch of the normal mode 
of the medium. 



3 THE PARAMETRIC INSTABILITY AND 
THE SOLITONS 

Assuming a uniform background field in the z-dircction and 
the variation of the field along this direction only, above set 
of Eqs. (4), (5) and (6) can be written as 

dp d{pv,) _ 

dt^ dz ~ ' 



dB 

dt 



d d 



dz 



1 d 



pdz S-TT 



d ^ d 

di^^'^d-z 



Bz dB 



An p dz 

i CiSh Po dB 
p dz 



(23) 



Here Ci — B-^/ ^/TiTpo, 5h ~ va/ojh is the neutral 
skin depth, C'l 5h = B^rin/B, B = B^ + iBy and 
V = Vx + ivy. A medium described by the above set of 
Eqs. (23) admits circularly polarized linear waves {B,v) = 
{Bo, Uq) expi (wo t—k^z), as an exact solution. This can 
easily be seen by linearizing Eqs. (23) around a static equilib- 
rium state. The linear dispersion relation provides a relation 
between the wave number fco and frequency wo of the wave 



2 

Wo ■ 



klCj 



1 ± 



LUH 



(24) 



The amplitudes Bo and (7o of the waves are related by 

-1 



Bo (jJo 
ko Bz 



1± 



(25) 



We note that the stability analysis of the steady-state back- 
ground consisting of the unperturbed as well as the circu- 
larly polarized pump waves, i.e. B = {Bx{z) , By{z) , Bz) and 
= (Ux{z), Uy(z), 0) with a constant density, indicates that the 
system is unstable, the so called modulational instability. 
The modulational instability has been cxtcusivcly studied in 
the past (HoUweg 1974; Stenflo 1976; Derby 1978; Goldstein 
1978; Brodin & Stenflo 1988; Pandey & Vladimirov 2007). 
Assuming such a steady-state background and linearizing 
Eqs. (23) one ends up with the following set of equations. 



dSp dSuz 

-m+'-d^ 



0, 



~9r 



■ ujo 5u- 



Bz fdSB^ 



An p \ dz 



+ ko SB. 



dSu- 



dt 



Bz 



'dSB- 



47rp \ dz 

dSuz 
"dT 



+ iJo Su+ — ko Uo Suz 



Bz Bo ; 
+ -y—T op , 



■ ko SB. 

gl dSp 

p dz in p 



dSB^ 
dt 



LOO 5B- = Bz (^^^^ + ko 5u- 



A-K p'^ 

Bo d5B+ 
dz ' 
dSuz 



B< 



dz 
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dz 



dSB- 
dt 



+ wo 5B+ = B, 



fdS-L 



V dz 



+riH 



d^5B+ 



+ 2 



tIh ko Bo dSp 
p dz 

ko Su+^ + ko Bo 5uz 

dSB- 



dz 



k'oSB+ , (26) 



While writing the above set of equations, we have assumed 
that SBz — 0. Here 5B+ — 5Bx cos(f> + 5By sintj), SB- = 
5Bx sin (f> — SBy cos <j) , 5u+ = Sux cos (f> + 5uy sin <j) , Su- = 
5ux sin 4>—duy cos(/>, and fiH = B^rin/B. Fourier analyzing 
the above set of equations (26) , we will end up with following 
8"* order dispersion relation. 



uj + aru! +... + oia; + ao = 0. 



(27) 



The coefficients ar, 
ing UJ = Lu/iJo, X 



a = die], k'ci 



Qo arc given in the Appendix. Defin- 
kjko, b = ijJo/luh, F± = 1/(1 ±6), 
= uJoX^ PF, a = Bo/Bz and noting 
that since the transition from stability to instability pro- 
ceeds through w = 0, in the w — » limit, we see that 
ao 



(28) 



In the long wavelength limit retaining only ~ 0{k) terms in 
the coefficients oi and oo one may write 



[F±i-(l + 6)] -2F±ix, 



(29) 
(30) 



Since for the left-circularly polarized pump waves, F+ = 
1/ (1 + b), from Eq. (30), ao = 0. Thus in view of Eq. (28), 
we shall anticipate that the instability should decreases in 
the vicinity of fc = 0. The numerical solution of the dis- 
persion relation (27) for the left circularly polarized (Fig. 
3(a)) indeed displays this feature. The normalized growth 
rate is plotted against the normalized wave number for 
uio/ijJH = 0.01,0.1 and 0.15. We see that the growth rate 
decreases in the vicinity of fe — » 0. With the increase in 
ujo /loh , the growth rate increases. As is clear from the wave 
dispersion relation (24), the increase in cjo/ujh implies the 
increasing importance of the Hall term (J x B), which ap- 
pears due to the relative drift between the plasma and the 
neutral (i.e. ion Hall-/?i ^ 1). This drift is caused entirely by 
the coUisional coupling of the ions with the neutrals. There- 
fore, the increase in uo/oJh implies the enhanced availabil- 
ity of the free energy reservoir. Hence with the increasing 
ujo/ijJh, growth rate increases. 

In Fig. 3(b) we plot the real (dotted line) and imaginary 
(solid line) part of the frequencies for luq/loh = 0.01. We see 
from Fig. 3(b) that the growth rate of instability causes 
the sharp decrease of the real frequency. It suggests that 
most of the available free energy of the pump is invested 
in the growth of the fluctuations resulting in the decay of 
the real part of the frequency. However, with the decrease 
in the growth rate, most of the the pump energy goes to 
the real part of frequency, as seen from the figure. Further, 
we note that before the real part of frequency shows signs 
of recovery, the imaginary part start to fluctuate and finally 
drops to zero. 

In Fig. 4 the growth rate for the left circularly polarized 
wave is plotted against the wave number, for different values 




S 



(b) |3 = 0.1 




Figure 3. The growth rate of the left circularly polarized 
wave Im[uj]/uJH against k/ko is shown in Fig. 3(a) for 
luq/ujh — 0.01 ,0.1 and 0.15 with plasma f3 = 0.1 and 
Bo/B^ = 0.1 . In Fig. 3(b) both real (dotted hne) and 
imaginary (solid hne) part of the frequency is given for 
'jJo/'^H = 0.01. All other parameters are same as in 
Fig. 3(a). 



of plasma beta, /3 = 0.01 , 1 , and 5. Although the growth rate 

is dependent upon the plasma compressibility, the depen- 
dence is not linear. For example, when plasma l3 increases 
from 0.01 to 5, the growth rate almost doubles. However, 
when both acoustic and Alfven speed equals, i.e. /3 = 1, the 
growth rate is highest. For (5 — 10, the growth rate is similar 
to the case when l3 — 5. Therefore, the highest growth rate 
is when va ~ Cs- The growth rate is not very sensitive to 
the the variation of the pump amplitude, Bq/Bz and hence 
is not given here. 
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Figure 4. Same as in Fig. 1, with varying plasma /3. 

For the right-circularly polarized waves, when F_ = 
1/(1 — ujo/loh), the waves can grow at a rate 

Im[uj] ^ F- , (31) 

in the vicinity of a; = 0. Since the growth rate (Eq. 31) 
is inversely proportional to (1 — ujo/loh), near ujq ~ loh, 
when wave frequency matches the Hall frequency, the in- 
stability can grow resonantly. Therefore, the growth of the 
instability is quite different for the left and right-circularly 
waves. Whereas for the left-circularly polarized waves, the 
instability decreases in the neighbourhood of = 0, for the 
right-circularly polarized pump, the instability may become 
large near fc = 0. However, for the right circularly polarized 
waves, unbounded growth of the instability is not possible 
since linear approximation will brealc down once fluctuation 
becomes comparable to the background equilibrium quanti- 
ties. 

In Fig. 5(a) the growth rate and in Fig. 5(b) correspond- 
ing real part of the frequency is shown for the right-circularly 
polarized waves. When loq/loh = 0.9, the growth rate be- 
comes very large. The free energy is resonantly pumped into 
the fluctuations with increasing ujo /ujh ■ The physical system 
behaves more like a driven oscillator. The resonant driving is 
indirectly related to the neutral-plasma collisions. The rela- 
tive drift between the plasma and the neutral causes a Hall 
field over the Hall time tnitH = ^^h^)- ^^'^ Alfvonwave 
propagation time uJq^ becomes comparable to the Hall time 
tu, the energy is freely fed to the fluctuation by the pump to 
the plasma particles. Resulting free energy causes the large 
growth rate. This behaviour can be seen from the analyti- 
cal expression Eq. (31). It should be noted from the corre- 
sponding curve in Fig. 5(b) that the real frequency shows a 
sharp decline for wq/uih = 0.9. This suggest that almost all 
the free pump energy is resonantly used in the fluctuation 
growth. Similar behaviour is also noted for wo/wir = 0.95. 
Since the instability growth rate in this case is larger than 
when ojq/ojh = 0.9, the part of the available free energy 




Figure 5. The dependence of the growth rate on the 
amplitude of the right handed circularly polarized pump 
wave is shown in Fig. 5(a) and corresponding real part 
of the frequency is shown in Fig. 5(b) for lu/loh = 0.9, 
and 0.95. The plasma /3 = 10 and Bq/B^ = 0.1. 

drains faster from the real part of the frequency. This is in 
agreement with the well known behaviour of the oscillators 
near resonance although present system is more complex. 

It is believed that the granulations or convective mo- 
tions are responsible for the Alfvonwave generation in the 
photosphere. However the medium is capable of exciting 
very low frequency fluctuations only, i.e. when the wave fre- 
quency is much smaller than the collision frequency. Then 
the neutrals are dragged along by collisions and as a re- 
sult, they participate in the collision (Tancnbaum & Mintzcr 
1962). However, excitation of the high-frequency ideal MHD 
Alfvenmode is very unlikely. The high-frequency ideal MHD 
mode will be damped by collisions of plasma particles with 
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the neutrals as well as collisions of plasma particles with each 
other (Tanenbaum & Mintzer 1962; Petrovic et al.; Vranjes 
et al. 2007, 2008). 

The excitation of very low frequency Alfvenwave is 
a promising candidate for heating and acceleration of the 
solar plasma from coronal holes. Taking loh = 1 Hz for 
Wei ~ 10* Hz we see that the growth rate of the left-circularly 
polarized wave 0.3 wo suggests that the instability of the 
Alfven wave could be relevant to the turbulent heating. Since 
iOH = 10 — 10''Hz, the resonance condition for the right cir- 
cularly polarized wave implies that ujo must be in the same 
range. Therefore, it is quite likely that the low frequency 
right circularly polarized waves may resonantly extract en- 
ergy form the ambient surrounding. 



4 DISCUSSION AND SUMMARY 

The solar photosphere is a weakly ionized gas with neutral 
hydrogen carrying the inertia of the fluid. We show that the 
non-ideal effects are of paramount importance for the ex- 
citation and propagation of the waves. The hydromagnotic 
waves in such a medium arc an outcome of the balance be- 
tween the inertia of the neutral component with the defor- 
mation of the magnetic field, to which neutrals axe tied only 
by the ion-neutral collisions. Therefore, very low frequency 
Alfven wave can bo excited in such a medium. The high 
frequency Alfven mode (in the ion fluid) will be damped in 
such a medium (Vranjes et al. 2007, 2008). 

The Hall diffusion dominates the ambipolar and Ohmic 
diffusion in photosphere. The Hall scale Lh = (p/Pi)^^^ Si, 
is typically two order of magnitude larger than the usual 
ion-inertial length 5i = VAij^d (Pandey & Wardle 2008). 
The presence of the Hall effect will cause the excitation of 
the ion-cyclotron and whistler waves in the medium. 

We show that in an inhomogeneous photosphere, circu- 
larly polarized whistlers, are an exact solution of the ensuing 
equations. These waves can act as a pump wave. The result- 
ing parametric instability can make the medium turbulent. 

We note that derivative nonlinear Schrodinger equation 
(DNLS) equations can be derived from Eqs. (23) using re- 
ductive perturbation technique (Kennel et al. 1988) 



0, (32) 



where &x = B/Bz and a = va Cij^ (C|C^). One can show 
that Eq. (32) admits magnetically compressive (bright) and 
magnetically depressive (dark) solitons (Hada et al. 1998). 

Defining 5i,2 = (feo =F \/277^) with e = [VIva - Ci) /Ci, 
describing the shift of the travelling wave speed V from the 
intermediate speed and ri = ^ — tT where = x — Cit and 
T = Cit, the solutions can be written as 

o o 1 S2 — So 

J — JO ~i ITT r , 

1 + Tb sinh^ {a Kr]) 

3> = 3arctan (Cs tanh (aKij)) + 

sgn ( — — — ) arctan (Cs tanh (a /fry)) , 
V 2 a J 



and C*B = ^/sUS^Cb- Here < Si < 5*0 < 5*2. 

So — Si 



S = So- 
$ = 



1 + Td sinh^ {aKrj)' 
-3arctan(CD teuih {a K rj)) 

—sgn — — ^ arctan (Co tanh (a 77)) , 



S'z - SI 
S2 — So 



Cb 



So — SI 

S2-S0' 



(34) 



and Ci, = ./S^JsICd with K = 0.5^{S2 - So) {So - Si). 
The value 61 = represents the minimum of the trans- 
verse field magnitude in a refractive, dark soliton and 62 = 
\/S2 is the peak value of b± in a compressive, bright soli- 
ton. The magnetic field depression in dark maximizes when 
61 = 0, i.e. when transverse components by and bz vanish. 
The magnitude of the field inside the dark soliton is deter- 
mined by Bz — Bo cos9, implying that a large propagation 
angle will cause large magnetic cavities. The electric field for 
the magnetic cavities can be estimated as 



(35) 



BoVA { boT \l ] /c sin 9 . 



Tb = 



SI 



So — Si 



Cb = 



SO 



So — Si 



(33) 



Here c is speed of light. Such a field is available for the 
charged particle acceleration. We speculate that if such a 
soliton induced field indeed exists in the solar photosphere, 
then with increasing height, since Alfven speed will increase 
from 5.4 x lO^cms"^ at ft. = in the photosphere to 
4 X lO^kms"-^ at /i = 10^ km for a lOG magnetic field, 
the electric field E will vary typically between 0.15 V/m to 
1.2 V/m. Such a field may not be sufficient for accelerating 
the particles. However, if the field strength is increased by 
two to three orders of magnitude (as could be the case in the 
sunspots and pores) a million Volt field can easily exist that 
will be available for the acceleration and subsequent heat- 
ing. We do not know of any observational support for the 
existence of soliton in the photosphere at this stage and thus 
the hypothesis of soliton generated field and ensuing charged 
particle acceleration should be treated with caution. 

The following are the itemized summary of the present 
work. 

(1) Non-ideal MHD effect, namely Hall diffusion is im- 
portant in the solar atmosphere. 

(2) The circularly polarized whistlers can be easily ex- 
cited in the medium. 

(3) The photosphere can become parametrically unsta- 
ble. 

(4) The nonlinear derivative Schrodinger equation de- 
scribes the finite amplitude fluctuations. We speculate that 
the dark and bright soliton solutions may cause the particle 
acceleration in the medium. 
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APPENDIX A: THE COEFFICIENTS OF THE 
DISPERSION RELATION EQ. (18) 

Defining b = ojo/ojh, a = B/Bz and x = fc/fco, the coeffi- 
cients of dispersion relation are given as 

08 = 1 , a7 = 3 fe F± X , 

-x^ {13 + -2b'^F±) , 
-^^2bF±{l+x') + {a^ + 2) [l-bF± [l - x^)] 

+l-a^ -3b-bF±x\l3 + a'^)-b [l + F± (l + x^)] 

+^i^-2b [F±(l + x') + l + F±x' {(^ + a)] -(Al) 

= (/? + a') x" + F± x' [B (1 + x') + a'] + x" 

+2bx [F±x (1 - a^) - bxF±-bx^ fI (/3 + a^)] 
+bxF± (l + x-2bx)-2x'^ F± {a^ + 2) 
+b (1 + x^) [1 - bF± (1 - x^) -F± {1 + x^)] 
+ [l + x^ + F± 1 +x^ {l3 + a')] [1 + F± (1 + x')] 
-[l+x'^ (l-a^) - 6(1 + a;^) + F^^ - b (3 F± x'^] 

X (1 + a;^) [1-6F± (1-a;^)] . (A2) 

^=x{a^ + 2) [-F± 1 + 6 (1 _ + (1 + a:^)] + 2 x 

[1 + x^ {1-a'') -b{l+x^) {l + l3F±x^)+F^^] 
- [F± ^ - 6 (1 - x^)] [F± x^ {213 + o^)) +2x0^] 
+x^ {l3 + a'^){b-l) + 2bx [x^ {d + a'^) + F± x^ 
X (1 + x^) + a" x" (F± + 1)] - X [F± ' + 1 + x^] 
X [1 - 6 - - bF±x^ {f3 + a)\ + \2bx- 1 - x\ 

X [1 + a;' + Fg^ + (/3 + a')\ . (A3) 

||=[F^^-6(l-x^)-l-x^] 

[l-Fx' (1 - o') - 6 (1 -h x^) (l-F/3F±x^) 
+2x'^ \F±x^ (2/3 + a^) +2a^] + {\-bF± (l - a;^)] 
X - b^x" (1 + x^) + x" {13 (1 + x'^) - a^)] 

+2bF±x^ (6-1) {13 + a) -x^ [Fg^ + l + x^] 
X [/3 + a' + F± {13 (1 + x") + a") + a'] 
-X [26x-l-a;] [l - a - b - bF±x'^ {P + a)] ■ (A4) 

^=[Fg'-b{l-x')-l-x']x 

[x^ (2/3 + a^) +2a^] -2F±x [Fg'^ -bl3x^ + x^) 
W {& (l + x^) + a^)] - x^ (6 - F± x) {l3 + a^) 
X [l + F±(l + x^)] +x^ {l+x-2bx) X 

[0 + 0" + F± {13 (1 + x^) + a ) + a'] . (A5) 

|£-=F±a;^ (6-1) {f3 + a^) {l + x~2bx) 

-[Fg^-bf3x^ {l + x'')+x^ {13 {l+x^)+a^)] 

X [1-6F± (l-x^) -F± (1 + a;^)] . (A6) 



